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ABSTRACT 

When  the  ocean  seabed  is  considered  to  be  rigid,  the  ideal 
waveguide  model  can  be  used  to  model  the  shallow  water  en¬ 
vironment.  However,  a  more  realistic  ocean  waveguide  model 
treats  the  ocean  floor  as  a  boundary  between  two  different 
fluid  media.  In  this  paper,  a  frequency-domain  characteri¬ 
zation  of  shallow  water  environments  is  proposed  based  on 
this  realistic  waveguide  model  with  a  fluid  boundary.  First, 
the  time-frequency  characteristics  of  this  model  are  studied 
as  well  as  the  impact  of  the  environment  parameters  on  the 
dispersive  phenomena.  Then  a  frequency-domain  matched 
filter  receiver  is  designed  to  obtain  time-dispersion  diversity 
once  we  separate  the  modes  using  warping  techniques  in  the 
time-frequency  plane.  Simulations  demonstrate  that  the  new 
receiver  design  improves  the  bit  error  rate  performance. 

Index  Terms —  Time-frequency  analysis,  dispersive  chan¬ 
nels,  time-frequency  mode  separation,  underwater  acoustic 
communications,  time-dispersion  diversity 

1.  INTRODUCTION 

As  a  waveform  is  transmitted  in  the  shallow  water  environ¬ 
ment,  it  can  interact  with  the  ocean  bottom  and  surface,  re¬ 
sulting  in  dense  dispersion  due  to  the  water  medium  and  time- 
varying  (TV)  changes  of  the  ocean  environment.  Specifically, 
the  shallow  water  acoustic  environment  is  a  linear  TV  dis¬ 
persive  system  that  shifts  different  frequencies  by  different 
amounts  in  time  [1]. 

The  dispersive  effect  introduced  by  shallow  water  envi¬ 
ronments  can  severely  limit  the  performance  of  underwater 
acoustic  applications  such  as  sonar  and  communications.  It 
can  be  shown  that  the  propagation  characteristics  of  an  en¬ 
vironment  can  provide  us  with  a  specific  nonlinear  function 
that  defines  the  dispersion  in  this  environment.  Thus,  the 
corresponding  signal  representation  based  on  a  normal  mode 
model  can  be  applied  to  exploit  the  potential  diversity  sug¬ 
gested  by  the  model  when  the  acoustic  receiver  is  appropri¬ 
ately  designed  to  match  this  nonlinear  function. 
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In  [1],  a  shallow  water  characterization  was  considered 
that  matched  dispersive  signal  transformations  on  the  trans¬ 
mitted  waveforms  and  was  used  successfully  for  shallow  wa¬ 
ter  communications  to  obtain  time-dispersion  diversity.  How¬ 
ever,  this  characterization  was  only  applicable  to  signals  with 
very  high  bandwidth  as  it  assumed  that  the  transmitted  wave¬ 
form  was  an  impulse. 

In  [2],  we  presented  a  generalized  characterization  that 
was  applicable  to  a  larger  class  of  signals  by  using  the  normal- 
model  in  [3]  assuming  perfect  waveguide  conditions  (homo¬ 
geneous  fluid  layer  with  a  soft  top  and  rigid  bottom).  In 
this  paper,  we  extend  our  characterization  when  these  condi¬ 
tions  are  not  perfect  by  following  the  waveguide  model  with 
pressure-release  surface  and  fluid  boundaries  in  [4].  Using  the 
new  model,  we  propose  a  new  time-frequency  receiver  design 
which  operates  on  individual  modes  of  the  multi-component 
received  signal  instead  of  on  the  whole  signal.  This  upgraded 
receiver  design  of  an  underwater  communication  system  is 
expected  to  achieve  better  bit  error  rate  (BER)  performance 
for  time-dispersion  diversity  than  the  one  obtained  in  [2] . 

The  paper  is  organized  as  follows.  In  Section  2,  we  for¬ 
mulate  the  Pekeris  model  for  shallow  water  environments.  In 
Section  3,  we  present  the  time-frequency  (TF)  characteristics 
of  this  model,  and  we  investigate  the  impact  of  the  environ¬ 
ment  parameters  on  the  dispersion  effect.  In  Section  4,  we 
use  a  TF  mode  separation  technique  to  separate  each  com¬ 
ponent  of  the  received  signal  and  then  design  a  new  receiver 
structure  with  a  corresponding  optimal  detector.  In  Section  5, 
numerical  results  of  BER  performance  illustrate  our  improved 
performance  in  terms  of  diversity  order. 

2.  MODELING  FOR  SHALLOW  WATER 
ENVIRONMENT:  PEKERIS  WAVEGUIDE 

The  normal-mode  model  treats  the  ocean  as  a  waveguide  with 
plane,  parallel  boundaries,  representing  the  acoustic  field  in 
the  ocean  medium  as  a  sum  of  normal  modes.  We  first  discuss 
the  normal- mode  model  for  shallow  water  environments  with 
pressure-release  surface  and  a  fluid  seabed  following  [4,  5,  3]. 
A  simple  waveguide  model  of  the  ocean  is  shown  in  Fig.  1 
using  the  coordinate  system  (r,z),  where  Medium  I,  II  and 
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Medium  I:  air 


Fig.  1.  Waveguide  model  with  point  source  in  Medium  II  at 
r  =  0,  z  =  zo,  D  m  deep.  This  corresponds  to  a  pressure 
release  surface  with  a  fluid  seabed  environment  [4] . 


Ill  correspond  to  air,  ocean  water  and  seabed,  respectively. 
An  omnidirectional  point  source  with  signal  spectrum  X(f) 
is  located  in  Medium  II  at  r  =  0  and  z  =  zq.  The  ocean  is 
D  m  deep.  We  consider  the  sound  speed  in  Medium  II  as  a 
constant  c  m/s  and  density  p  kg/m3,  and  in  Medium  III  as  a 
constant  cb  m/s  and  density  pB  kg/m3. 

The  ocean  surface  (at  z  =  0)  is  modeled  realistically  as 
an  ideal  pressure  release  boundary  and  the  ocean  bottom  (at 
z  =  D)  is  modeled  as  a  boundary  between  two  different  fluid 
media.  The  normal  mode  model  is  given  by  the  solution  of 
this  ocean  waveguide  problem,  which  is  determined  by  the 
environment  parameters  and  satisfies  all  boundary  conditions, 
including  the  boundary  condition  at  the  source. 

After  a  detailed  derivation  in  [4]  and  [5],  the  received  sig¬ 
nal  spectrum  excited  by  X(f)  at  location  (r,z)  is  given  by 
the  Pekeris  waveguide  model: 


Np- 1 

Ypekeris{f)  =  X(f)  £  C„(/)0„(/).  (1) 

71=0 

Without  the  assumption  of  ideal  waveguide  condition  as  in 
[2],  the  nth  mode  is  characterized  by 
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is  the  wave  number  of  the  nth  mode  and  Np  is  the  largest 
mode  number.  The  parameter 


Cn{f)  =  A2n(f)sm(kzn(f)z0)sm(kzn(f)z) 


is  a  function  of  frequency  where 
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3.  TIME-FREQUENCY  CHARACTERISTICS  OF 
SHALLOW  WATER  ENVIRONMENT  MODEL 

The  TF  characteristic  of  the  acoustic  signal  is  determined  by 
the  modal  group  velocity  (MGV)  gn(f ),  due  to  the  propagat¬ 
ing  delay  of  the  frequency  /o  in  the  nth  mode  is  determined 
by  Tn(/o)  =  y~Uo)'  ^he  m°dal  SrouP  velocity  is  shown  in 
Fig.  2.  From  Fig.  2,  we  notice  that,  for  the  waveguide  model 
with  fluid  seabed,  the  MGV  approaches  c  when  the  frequency 
approaches  infinity,  and  approaches  cb  when  the  frequency 
approaches  the  cutoff  frequency  of  this  mode. 


Fig.  2.  Modal  group  velocity  (MGV)  function. 

We  analyzed  further  the  TF  characteristics  of  the  received 
signal  using  the  spectrogram  TF  representation  (TFR)  and  in¬ 
vestigated  the  effect  of  changing  the  parameters  of  the  envi¬ 
ronment.  Fig.  3(a)  represents  the  spectrogram  for  three  modes 
with  environment  parameters:  D  =  100  m  and  r  =  15  km. 
Fig.  3(b)  uses  the  same  set  of  parameters  except  that  r  =  30 
km.  Fig.  3  shows  each  mode  as  a  dispersive  curve  in  the  TF 
plane.  We  note  that  the  curves  stop  at  the  cutoff  frequency 
with  no  asymptotical  behavior  to  the  cutoff  frequency.  Fig.  3 
shows  that  longer  range  between  the  receiver  and  transmitter 
causes  more  dispersion  in  the  received  signal.  Fig.  4  shows 
that  with  larger  ocean  depth,  the  modal  cutoff  frequency  is 
lower,  and  with  a  smaller  depth,  the  dispersive  effect  is  more 
visible  at  higher  frequencies. 


Fig.  3.  Spectrograms  with  environment  parameters  (a)  r  =  15 
km,  D  =  100  m,  and  (b)  r  =  30  km,  D  =  100  m. 


Fig.  4.  Spectrograms  with  environment  parameters  (a)  r  =  15 
km,  D  =  30  m,  and  (b)  r  =  15  km,  D  =  100  m. 


4.  SYSTEM  DESIGN  FOR  COMMUNICATION 
DIVERSITY 

According  to  the  Pekeris  waveguide  in  (1),  the  received  sig¬ 
nal  spectrum  is  the  superposition  of  multiple  normal  modes. 
Considering  each  mode  as  a  communication  subchannel,  the 
time-dispersive  channel  can  be  interpreted  as  a  superposition 
of  subchannels,  and  using  their  TF  characteristics,  we  can 
separate  each  subchannel  and  design  the  corresponding  re¬ 
ceiver  to  exploit  time-dispersion  diversity. 

4.1.  Separation  of  normal  modes 

Fig.  3  illustrates  that  each  mode  of  the  received  signal  appears 
as  a  dispersive  curve  in  the  TF  plane,  which  makes  the  sepa¬ 
ration  of  the  modes  possible.  In  this  paper,  we  use  a  two-step 
TF  mode  separation  technique  based  on  warpings  to  separate 
the  modes.  The  two  steps  are  described  as  follows. 

During  the  first  step,  we  design  Np  —  1  instantaneous  fre¬ 
quency  (IF)  curve  separators  in  the  TF  plane,  represented  as 
en(£),  n  =  0,  •  •  •  ,  Np  —  2.  An  IF  curve  separator  is  a  curve 
situated  between  two  successive  modes  in  the  TF  plane  as 
shown  in  Fig.  5(a)  for  the  example  in  Fig.  3.  In  this  prob¬ 
lem,  all  environmental  parameters  are  assumed  to  be  known, 
so  that  the  TF  structures  of  the  received  signal  are  also  known 


Fig.  5.  (a)  Spectrogram  of  received  signal  with  the  IF  curve 
separators,  (b)  separated  TF  modes. 


using  the  model  in  (1).  Knowing  the  TF  structures,  we  ap¬ 
propriately  set  M  TF  points  (■ U,  fa),  i  =  1,  2, . . . ,  M  in  the 
middle  of  the  space  between  two  successive  modes.  The  M 
points  between  the  nth  mode  and  the  (n  +  l)th  mode  will 
constitute  the  IF  curve  separator  en(t).  In  the  next  step,  the 
IF  curves  will  be  processed  with  the  received  signal.  Thus,  the 
separation  lines  must  have  the  same  duration  as  the  received 
signal.  To  ensure  that,  all  the  IF  curves  begin  at  the  starting 
point  of  the  TF  structure  and  are  extended  with  a  straight  hor¬ 
izontal  line  until  the  end  of  the  TF  window. 

During  the  second  step,  we  process  the  different  modes  in 
the  TF  plane  using  a  warping  technique  [6, 7]  based  on  the  use 
of  the  local  harmonic  convolution  operator  [6].  Assume  that 
the  received  noisy  signal  in  (1)  is  r(t)  with  Fourier  transform 
(FT)  R(f)  =  YPekeris(f)  +  W(f),  where  W(f)  is  additive 
white  Gaussian  noise.  We  then  obtain  the  warping  functions 
from  the  IF  curve  separators  as  Cn(t)  =  j~  en(u)du, 
where  fc>  0,  and  compute  the  corresponding  generalized 
FT  of  r(t)  as  [8] 

Mr(X)  =  (  r{i)d^e-^x^tUt,  (3) 
J  p  dt 

where  p  is  the  domain  of  time  of  the  warping  function 
and  A  is  real  and  unitless.  To  obtain  the  first  mode,  we  com¬ 
pute  the  inverse  generalized  FT  over  the  range  (-00  ,fc/fr] 
as  follows: 

rfc/fr 

r0(t)=  /  Mr{X)ej2vX<:oit)dX,  (4) 

7—00 

where  fa  is  a  normalization  frequency.  Then  we  subtract  the 
fist  mode  from  the  received  signal,  i.e,  r(t)  =  r(t)  —  r0(t), 
and  the  remaining  modes  are  contained  in  r(t). 


To  obtain  the  seconde  mode,  we  apply  the  above  proce¬ 
dure  to  f(t)  using  the  IF  curve  separator  e\  (t).  Repeating 
the  above  procedure  for  each  IF  curve  separator,  we  can  then 
separate  each  mode  from  the  received  signal. 

For  the  example  in  Fig.  5,  we  use  the  mode  separation 
technique  on  the  received  signal  excited  by  the  waveform 
X(f)  =  v7  f°r  Np  =  3  modes.  The  IF  curve  separators  are 
shown  in  Fig.  5(a)  and  the  separated  components  are  shown 
in  Fig.  5(b). 

4.2.  Receiver  structure 

Under  the  Pekeris  waveguide  assumption,  the  model  in  (1) 
can  be  seen  as  an  inherent  frequency  domain  transfer  func¬ 
tion.  However,  due  to  the  distortion  and  signal  propagation 
uncertainty  in  realistic  shallow  water  environments,  a  more 
suitable  model  of  random  signal  embedded  in  random  noise 
consists  of 

Np  —  l 

R(f)=Y(f )  +W(f)=X(f)  J2DnCn(mn(f)  +  W(f), 

n= 0 

(5) 

where  R(f)  is  the  spectrum  of  the  received  signal,  Dn  is  a 
random  process  incorporating  amplitude  and  phase  variability 
of  the  source,  and  W(f)  is  additive  white  Gaussian  noise  with 
variance  cr ^ . 

We  propose  a  filter  bank  receiver  scheme  based  on  this 
signal  model  which  can  exploit  the  frequency  domain  disper¬ 
sion  diversity  existing  in  the  normal  modes.  The  received 
signal  spectrum  R(f)  in  (5)  is  first  processed  as  described 
in  Section  4.1  to  separate  each  mode,  and  then  the  outputs 
are  combined  in  the  minimum-probability-of-error  sense  to 
obtain  the  transmitted  information  symbol  b.  This  receiver 
structure  is  illustrated  in  Fig.  6. 

We  design  the  transmission  waveform  as  X(f)  =  v7  f°r 
/  >  /o  where  /o  is  the  cutoff  frequency  of  the  first  mode,  and 
then  we  employ  the  corresponding  frequency  domain  matched 
filter  to  each  separated  mode  component  of  the  received  sig¬ 
nal.  The  nth  mode  component  obtained  from  the  separation 
processing  is  expressed  as 

Rn{f)  *  bX(f)DnCn(f)Qn(f)  +  Wn(f),  (6) 

where  Wn(f )  is  non- white  Gaussian  noise,  with  variance  Oyyn 
that  resulted  from  W (/)  when  the  nth  mode  is  separated. 

We  first  define  the  receiver  matched  waveform  Un(f )  = 
X(f)Cn(f)Qn(f)  for  the  nth  mode.  Then,  as  each  mode  is 
separated,  we  can  perform  individual  matched  filtering  on  the 
noth  mode  by  DnoUno(f );  this  can  avoid  interference  with 
the  other  modes.  The  output  of  the  noth  matched  filter  can  be 
expressed  as 

Zno  =  (Rno(f),DnoUno(f))  =  I  Rno(f)D*noK0(f)df. 

(7) 


R(f) 


Estimated  Symbol 


Fig.  6.  Receiver  structure  scheme. 

Using  the  expression  of  Rn(f)  in  (6),  we  can  rewrite  (7)  as 

Zno  *  b\Dno \2Juno (f)U*0 (f)df+J\Vno {f)D*noU*0 (f)df. 

As  the  output  of  this  filter  bank  is  given  by  the  concate¬ 
nation  [Zq  ...  ZNp_ i] ,  the  combination  in  the  minimum- 
probability-of-error  sense  is  defined  as  the  arithmetical  mean 
of  the  filtering  outputs.  Thus,  the  estimation  of  the  transmit¬ 
ted  information  symbol  can  be  expressed  as 

fNP-i  \ 

b  =  5R  I  I  '  ® 

where  Sft(-)  is  the  operator  that  takes  the  real  part.  Assuming 
an  antipodal  symbol  transmission,  i.e,  b  =  +1  orb  =  —  1,  and 
that  the  probabilities  of  transmitting  +1  and  —1  are  equal,  the 
minimum  error  probability  detector  can  be  expressed  as: 

decide  b  =  +l  if  b  >  0, 

decide  b=—  1  if  b  <  0. 

5.  PERFORMANCE  ANALYSIS 

As  the  modes  are  separated,  each  mode  can  be  considered 
as  a  subchannel,  and  the  aforementioned  matched  filter  re¬ 
ceiver  can  avoid  cross-correlation  interferences  between  dif¬ 
ferent  modes  of  the  received  signal.  Hence,  this  matched  filter 
receiver  can  exploit  the  potential  diversity  existing  among  the 
different  modes.  In  this  section,  we  first  derive  the  theoretical 
probability  of  error  for  the  aforementioned  receiver,  and  then 
we  obtain  the  BER  results  from  numerical  simulations.  We 
also  compare  the  performance  of  this  receiver  with  the  per¬ 
formance  of  a  receiver  that  processes  the  complete  received 
signal. 


Let  us  consider  the  matrix  of  frequency  domain  mode 
autocorrelations  defined  as  A  =  diag(A0:  A1: . . . ,  A/vP-i) 
where 

Ai  =  J  ui(f)u;(f)4f,  t  =  0, . . . ,  NP  —  1 

and  the  matrix  of  modal  noise  variance  a  =  diag(a^v.)  for 

i  =  0,  ...,NP  —  1.  Defining  ^  =  cr_1V/AQv//A^  and 
m  =  \4Fita  e[d],  where  E[  ]  is  the  expectation  opera¬ 
tor  and  Q  is  the  covariance  matrix  of  the  coefficients  vector 
D  =  [Do,  Di, . . . ,  DiVp-i]r,  we  can  derive  the  theoretical 
BER  as  the  probability  of  error  given  by  [2] : 


Perr  — 


/f 

r  y  1 

det  ^  + 1 

/  0 

sin2  0 

1  -1 


^—r\A  (y^+sin2  6  I)  1m 


(9) 

where  we  denote  T  and  f ,  respectively,  as  the  transpose  and 
conjugate  transpose  operations.  The  BER  obtained  from  the 
simulation  and  the  theoretical  expression  are  shown  in  Fig.  7, 
using  the  environment  parameters  D  =  100  m,  r  =  15  km, 
c  -  1500  m/s,  cb  =  1800  m/s,  p  =  1000  kg/m3  and  pB  = 
1800  kg/m3  (as  in  Figures  3  and  5).  The  results  show  the 
BER  as  a  function  of  the  signal-to-noise  ratio  (SNR)  for  the 
transmission  frequency  band  5  —  250  Hz. 

We  notice  from  Fig.  7  that  the  simulation  results  are  worse 
than  the  theoretical  results.  This  is  because  the  mode  separa¬ 
tion  technique  caused  loss  of  information,  especially  at  higher 
frequencies,  which  leads  to  less  accurate  matched  filtering. 
However,  if  we  assume  no  information  is  lost  in  the  proce¬ 
dure  of  separation,  i.e,  the  separation  can  generate  the  exact 
modes,  then  the  numerical  simulation  results  are  close  to  the 
theoretical  results. 

We  also  compared  the  BER  performance  of  the  new  re¬ 
ceiver  with  the  performance  obtained  in  [2]  where  the  matched 
filtering  was  applied  to  the  whole  multi-component  received 
signal.  Fig.  7  shows  that  the  new  receiver  performed  better 
according  to  the  numerical  simulation  results. 


6.  CONCLUSION 


Fig.  7.  Theoretical  BER  performance  of  the  new  receiver  (di¬ 
amonds),  BER  performance  of  the  new  receiver  obtained  by 
numerical  simulation  (squares),  BER  performance  of  the  new 
receiver  obtained  by  numerical  simulation  using  exact  modes 
(line)  and  BER  performance  of  the  former  receiver  (circles). 
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2066-2076,  May  2007. 

[2]  J.  Zhang  and  A.  Papandreou-Suppappola,  “Time- 
frequency  based  waveform  and  receiver  design  for  shal¬ 
low  water  communications,”  in  Proc.  IEEE  Int.  Conf. 
on  Acoust.,  Speech,  Signal  Processing ,  vol.  3,  pp.  1149- 
1152,  April  2007. 
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Space-Time  Signal  Processing ,  CRC  Press,  1999. 
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We  developed  a  frequency  domain  characterization  of  real¬ 
istic  shallow  water  environments  and  analyzed  its  dispersive 
characteristics.  Following  this  model,  we  developed  the  ap¬ 
propriate  receiver  design  based  on  using  a  warping  technique 
for  modal  separation.  This  receiver  structure  aims  to  exploit 
the  diversity  existing  in  the  system.  Simulation  results  show 
that  BER  performance  and  diversity  order  were  improved  by 
the  aforementioned  system  design. 
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